Zn-impurity effect and interplay of s±- and s ++ -pairings in Fe-based superconductors 
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We report theoretical and experimental studies of the effect of Zn-impurity in Fe-based supercon- 
ductors. Zn-impurity is expected to severely suppress sign reversed s± wave pairing. The experimen- 
tally observed suppression of T c under Zn-doping strongly depends on the materials and the charge 
carrier contents, which suggests competition of s++ and s± pairings in Fe-base superconductors. 
We study a model incorporating both s++ and s± pairing couplings by using Bogoliubov de-Gennes 
equation, and show that the Zn-impurity strongly suppresses s± pairing and may induce a transition 
from s± to s++-wave. Our theory is consistent with various experiments on the impurity effect. We 
present new experimental data on the Zn-doping SmFei-^Zn^AsOo.gFo.i of T c = 50K, in further 
support of our proposal. 

PACS numbers: 



I. INTRODUCTION 

One of the most important issues in the high T c Fe- 
based superconductors (FeSC) is their pairing symme- 
try Theories based on antiferromagnetic (AF) 
spin fluctuations have predicted s± pairing, where the 
superconducting (SC) order parameters on the hole and 
electron Fermi pockets have opposite signs [3, The 
proposed symmetry is consistent with a number of ex- 
periments, such as the spin resonance peak in neutron 
scattering Q, sensitive SC junction data plUf, and quasi- 
particle interference in tunneling experiments 0]. How- 
ever, the pairing symmetry in FeSC may not be universal, 
and there are evidences for different pairing structures as 
discussed in a recent review [1(| ■ 

The effect of disorder to the superconductivity is an 
important test to the pairing symmetry According to 
Anderson's theorem, the conventional s-wave supercon- 
ductivity is insensitive to non-magnetic impurities. The 
sign reversed s± superconductivity is, however, sensitive 
to non-magnetic impurities which scatter inter-band elec- 
trons. Replacement of part of Fe-atoms by Co or Ni 
in a parent compound of FeSC leads to superconductiv- 
ity. However, the role of the Co or Ni doping is more 
subtle and remains controversial. One scenario is that 
the doping introduces additional electron carriers. This 
scenario is supported by the angle resolved photoemis- 
sion spectroscopy, which indicates the shrinking of the 
hole pockets On the other hand, recent resonant 

photoemission spectroscopy and density functional cal- 
culations indicate that Co doping is covalent and intro- 
duces disorder [12]. ft is plausible that the Co doping 
introduces both carriers and disorder 131. Zn-ion has 



a 3d 10 configuration, hence a very high electric poten- 
tial to charge carriers. Replacing a Fe-atom by Zn in 
FeSC introduces inter-band scattering and is expected 
to severely suppress the s± superconductivity. There- 
fore, the Zn-doping is an effective test to the s± pair- 
ing in FeSC. There have been several experiments on 
the Zn-doping effect on FeSC, including so-called 1111 
compounds LaFei-^Zn^AsOi-yFj, [3, EEj], and more 
recently 122 compounds BaFe2(i- a; -; / )Zn2 2 ;Co2j / As2 and 
SrFei.8-2£cZn2:rCoo.2As2 pUf ■ The results are mixed at 
present, which appears to be strongly dependent of ma- 
terial and charge carrier content. The experimental data 
on the 1111 compounds may be divided into two cate- 
gories. The optimally doped LaFeAsOo.gFo.i [3 i s m ~ 
sensitive, but the over-doped LaFeAsO0.s5F0.15 is very 
sensitive to the Zn- impurities [l5j]. The effect of Zn- 
doping on Co-doped 122 compounds clearly shows the 
suppression of superconducting transition temperature 
T c , but the reduction is much slower than the theory 
predicted A careful examination indicates that the 
suppression of T c may be saturating at large Zn-doping 
to some of the compounds. Note that it is not easy to 
dope Zn into the Fe lattices uniformly even under high 
pressure, and reliable data is only available up to 6% Zn- 
doping at present. Therefore the experimental data are 
not complete. Nevertheless, the available experiments 
on Zn-doping indicate complexity of the effect, and sug- 
gest possible competition of sign changed s± and sign 
unchanged s++ pairings in FeSC. 

In this paper, we use a two-orbital model for FeSC 
including both on-site (or s++) pairing coupling go and 
next nearest neighbor (NNN) intersite (or s±) pairing 
coupling #2 to study Zn-impurity effect, which may help 
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understand the complex result of the Zn doping effect 
on 1111 and 122 compounds. We apply Bogliubov de- 
Gennes (BdG) equation to study the model on a finite- 
size system. The two SC pairings in the multi-band sys- 
tem show interesting interplay. They may mix but also 
compete with each other. The disorder strongly sup- 
presses the intersite pairing, and its effect to the super- 
conductivity depends on the strength of go- For large 
go, <72 plays little role and the pairing is s++ and is ro- 
bust against the disorder. For small go, the pairing is s± 
and the disorder strongly suppresses superconductivity. 
For moderate value of go, the disorder may enhance the 
on-site pairing and induce a transition from s± to s ++ su- 
perconductivity. We further study the interplay between 
go and gi in a clean system and show that the disorder 
effect on the gap functions is similar to the reduction of 
<?2- Our theory is consistent with the Zn-doped impurity 
experiments on 1111 and 122 compounds, and suggests 
multi-pairing couplings in some of the FeSC. We present 
our new experimental data of the Zn-impurity effect on 
the very high T c = 50K Sm-1111 compound. The lat- 
tice constant measurement show that the Zn- atoms are 
doped into the Fe-lattice uniformly up to 6%. The re- 
sults appear to indicate possible saturation of T c under 
the Zn doping, consistent with the present theory. 



II. MODEL HAMILTONIAN AND MEAN FIELD 
THEORY 




FIG. 1: (Color online) Hole (red) and electron (blue) Fermi 
pockets obtained in the two-orbital model Eq. (2). Points A 
and B are the representative k points for the hole and electron 
pockets, respectively. 



the Zn-impurity site, 



Hi mp — i ^2 e c\ a Cj a , 

2£imp a 



(3) 



where i sums over all the impurity sites, and we consider 
the large I case (7 = 24ii in the actual calculation) to 
model the large repulsion to an electron at the Zn site. 
The pairing Hamiltonian is modeled by 

#pair = - (^J C Lt C iW C ^^ Clc "t + h.C.), (4) 

<ij> 



We consider a model Hamiltonian 



H = H a +H 



pair 



H 



imp, 



(1) 



which includes a tight-binding kinetic term Ho, a pairing 
interaction i? pa ir, and a disordered term H lmp . For Ho, 
we consider a two-orbital model with d xz (orbital 1) and 



d yz (orbital 2) as proposed by Raghu et al. (l7j . 



Ho 



E 



(2) 



where C\ a = (cl 



and h 



t?4 is the electron 



hopping term between orbital a at site i and orbital /3 
at site j on a 2-dim. square lattice of Fe-atoms (lattice 
constant a = 1). While this model may be an over- 
simplified one to describe many detailed material prop- 
erties of FeSC, it should capture the basic feature of the 
disorder effect to the pairing in systems with multi-Fermi 
surfaces. The non- vanishing hopping matrix elements 



are if? = -/i, * 



li 



t 12 - - 

t,t-\-x-\-y 



t 21 - - = 

i,t-\-x-\-y 



= t. 



4-a.a j 

i,i+x+y ~ l 3, 



— ^2) ti,i+y — Hi 

-<4. We choose t\ = 1 as the en- 
ergy unit, and ^ = 1.6, t% = —1.3, £3 = t± = 0.85, which 
gives Fermi surfaces with hole pockets near the T- and M 
points, and electron pockets near the X- and Y points in 
an extended Brillouin zone as plotted in Fig. 1. 

We consider randomly distributed impurities on the 
lattice and introduce an on-site repulsive potential on 



where the pairing coupling Vij includes an on-site term 
go > and an NNN intersite term 172, 



(5) 



with t the vector of the two NNN site displacement. Note 
that go term favors s++ and 32 term favors s± symmetry. 

We introduce a mean field gap function Af" = 
Vij (cjaiCiai). Our calculations show that the inter- 
orbital pairing A] 2 is very tiny and will be neglected 
below. The BdG equation for the mean field Hamilto- 
nian then reads 



E 



foij ^ij 
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= E n 
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with A, 



«il,<7 
Ui2,cr 

tion is 



AijI, and / an identity matrix. u 4 ,, 



The self-consistent equation for the gap func- 
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For the form of V v] in Eq. ©, we define A% a [i) = A 
and Ar« = E ? A?^/4. 



a a 
ii i 
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III. NUMERICAL RESULTS 

We now discuss the numerical solutions of H. In our 
calculations, for each impurity content, the impurity po- 
sitions are randomly distributed and the statistical aver- 
ages are taken over 400 times. We consider three typical 
cases: (i) go is large and dominant; (ii) g 2 is large and 
go is weak; and (iii) gi is dominant but go is moderately 
large. In case (i), the SC pairing is always s ++ and the 
superconductivity is robust against the impurity as we 
expect from the Anderson theorem. 

In Fig. 2 (a) and (b), we show the spatially averaged 
on-site and NNN inter-site pairing amplitudes Ao = Ag a 
and A2 = A 2 Q as functions of the impurity concentra- 
tion nj mp for cases (iii) and (ii). Also shown are the gaps 
at the hole pocket (point A [(0, 0.227r)]) and at the elec- 
tron pocket (point B [(0.627T, 0)]), which are the Fourier 
transform of the impurity averaged gaps in real space. 
In the case (ii) of weak on-site pairing, the impurities 
strongly suppress A2 as shown in the Fig. 2(b). Ao is 
tiny and the SC gap functions A a and Ab monotoni- 
cally decrease as rii mp increases. Because of the finite 
lattice size, our study is limited to the short coherence 
length or the strong pairing coupling cases, which require 
?^imp ~ 0.15 to destroy the superconductivity. We expect 
this value to be much smaller in weaker pairing coupling 
cases. 

Case (iii) is most interesting, and our theory shows an 
impurity driven phase transition from s± to s++ pairings. 
In the absence of impurity, g 2 dominates and the pairing 
is s±. As shown in Fig. 2(a), the pairing symmetry 
remains to be s± at n lmp < 0.02, and the gap amplitudes 
on k points A and B are monotonically suppressed as 
n imp increases. At 0.02 < n imp < 0.05, |A 2 | decreases, 
and |A | increases. At n imp > 0.05, both A^ and Ab 
are positive and we have s++ pairing. It is interesting 
to note that the on-site pairing may be enhanced by the 
impurities due to the suppression of the NNN pairing. 

In the case of weak on-site pairing, as shown in Fig. 
2(b), the impurities strongly suppresses A 2 . We have ex- 
amined the SC order parameters in real space and found 
that the disorder does not result in severe pair-breaking 
effect to the on-site pairing measured by Aff, whose 
peak amplitude is almost unaltered by the impurities. On 
the other hand, the non-magnetic impurities not only de- 
stroy NNN SC pairing order parameter A"" + - + - in larger 
spatial areas, but also weaken the peak amplitude of the 
SC pairing immensely. 

As we have demonstrated, the impurities suppress the 
NNN pairing order parameter A 2 . This effect is similar 
to the reduction of gi in the clean sample. Therefore, 
tuning n lmp in the disorder system is similar to tuning 
gi in a clean system [ijj]. Below we shall study SC or- 
der parameters and T c in the model Hamiltonian H as 
functions of g 2 in the absence of disorder to mimic the 
impurity effect. This enables us to further reveal the 
interplay between the SC pairings of s++ and s±. 
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FIG. 2: (Color online) Upper panel: The gap functions at 
hole and electron Fermi pockets Ah and A e as functions of 
impurity density ni mp , obtained in the mean field solution for 
H. (a): for modestly strong on-site pairing coupling go = 1.8; 
(b): for weak on-site coupling go — 0.8. Insets: spatially 
averaged gap functions Ao (on-site) and A2 (NNN inter site). 
In both cases, the NNN coupling g-2 = 1.6. Lower panel: 
SC gaps calculated by the simplified BCS formalism, with 
ATe(O) = 0.12, N h (0) = 0.1, and tv D = 0.8. (c): g = 1.8; (d): 
go = 0.8. 



IV. T c REDUCTION: THEORY AND 
EXPERIMENTS 

For a clean system, we have lattice translational sym- 
metry, and the gap function Eq. (7) becomes 



A m(k) = -^y m „(k,kO tanh (f" k//2) A n (kO 



(8) 



where 



m, n 



are the band indices, -E„k = 
e„(k) 2 , e„(k) is the single particle en- 



VA^F 

ergy. The summation is taken only in the vicinity of 
Fermi pockets with an energy cut-off u>d- The pairing 
potential V m , n (k, k') describes the coupling between gap 
function on various Fermi pockets, and with Eq. (U), we 
have 

K™(k,k') = J2u ma (-k)U ma (k)U na (k / )U na {-k l ) 

OL 

x (50 + 4g 2 cos q x cosq y ), (9) 

where m, n are band indices, a is orbital index, U (k) is 
the transformation matrix between bands and orbitals. 

In our two-orbital model, there are four Fermi pockets, 
two for hole bands at T and M points respectively and 
two for electron bands at X and Y points respectively, 
which makes it very difficult to solve Eq. (0) analytically. 
So in the following, we will ignore the size of the pockets 
and assume there are four point-like Fermi surfaces at 
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r, X, Y and M with finite density of states. And we 
also assume the summation in Eq. (|5J) are only over the 
four momentum T = (0,0), Y = {ir, 0), X — (0, 7r), and 
M = (tt, vr). 

Then we consider the transformation matrix under this 
approximation. In the two-orbital model, the two or- 
bitals, d xz and d yz , mixes strongly in the hole Fermi 
pockets. On the other hand, the two orbitals can 
be connected by a C4 rotation. So in the case of 
point-like hole Fermi surface, it is obviously that the 
two orbitals contribute equally to the hole pockets, i.e. 
Uh,xz(yz)[F(M)) = where h denotes the hole band 
and xz and yz denote the two orbitals. On the other 
hand, the two electron pockets are dominated by d xz 
and d yz orbital respectively. So under the small pocket 
approximation, we have U eiXZ (Y) = U e ^ yz (X) = 1, and 
U e y Z (Y) = U exz (X) = 0. And the nonzero pairing po- 
tentials are V^(r,T) = V hh (T,M) = V hh (M,M) = f , 
V ee (X,X) = V ee (Y,Y) = wo, and V he [T(M),X{Y)} = f , 
where v = go + 4g 2 and v 2 — go ~ 4g 2 . 

In the small pocket approximation, the gaps on the 
two electron pockets should be same because of the C4 
rotational invariance of the iron pnictide. Though the 
gaps on the hole pockets may be different, we still assume 
they are equal for simplicity. So with the above pairing 
potentials, we can solve the gap equation ([5]) and get the 
critical temperature 



with 



k B T c = l.lAuj D e 



-1/N h (0)v 



(10) 



A) (So + 
where 



y/(l + A) 2 (<7o + 4<? 2 ) 2 - 16A.g <? 2 ], where A 
N e {0)/N h (0), and A e (0) and N h (0) denote the density 
of state at the Fermi level of electron and hole pockets 
respectively. 
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FIG. 3: (Color online) Left: The critical temperature T c as a 
function of <?2 with strong, weak, and moderate on-site pair- 
ing coupling go- Right: Zn- impurity effect on T c in various 
Fe-based superconductors observed in experiments. The ref- 
erences of the data are listed. 

Our calculation on the critical temperature for vari- 
ous impurity concentrations, depicted at the left panel of 



Fig. [3l reveals that the different impurity-doping behav- 
iors observed in FeSC [HI, [l6| may be characterized by 
the strength of the effective on-site pairing potential go- 
There are three types of cases for the disorder effect. In 
the case of large go, where the on-site pairing dominates, 
T c is hardly suppressed by the Zn doping. In the case of 
weak go, superconductivity is destroyed by the impurity. 
When go is comparable with g 2l as Zn impurity concen- 
tration increases, T c is initially suppressed rapidly and 
then saturate. The experimental facts seem to support 
the above scenarios and the effect of Zn doping depends 
on the material and the charge carrier concentration . In 
LaFei-zZn.AsOo.gFo.i (Ref. 0), T c are insensitive to 
the Zn- impurity, and may be explained due to large go . In 
the over-doped LaFej^Zn^AsOo.ssFo.is (Ref. [lj|) and 

BaFe 



LaFeAsOo.85 (Ref. [19(), in BaFe 2 (i_ a: _ y )Zn2 a; Co23 / As2 
(Ref. ;i6j), and in LaFei-^CoyZn^O (Ref. 0), T c 
decreases rapidly with the Zn doping, and may belong 
to the category of weak go- In SrFei.g-2a;Zn2a:Coo.2As2 
(Ref. 16]). T c was found to decrease slowly and has the 
tendency to saturate although higher Zn-doping will be 
needed to confirm the speculation. These scenarios are 
summarized in Fig. [31 which shows the critical tempera- 
ture vs gi at different go compared with the experimental 
data of the three types of materials that behaves differ- 
ently upon Zn doping. 

The moderate value of g 2 case is most interesting, 
for it reflects the competition between the two SC pair- 
ings. To further explore this possibility, we have pre- 
pared SmFei-^ZnzAsOo.gFo.i system with Tc = 50K 
and studied systematically the Zn-impurity effect to T c 
experimentally. The results are summarized in Fig. 4. 
We have measured the change of the lattice constant due 
to Zn-doping and confirmed that Zn-atoms are indeed 
doped into the iron sites up to 6% of Zn doping, see Fig. 
4(b) [2l|. Beyond this doping, our data indicate that 
some Zn-impurity may not enter into Fe-lattice so the 
measurement of T c may not correspond to the uniformly 
doped Zn-impurities. The main experimental result of T c 
versus Zn concentration on this very high T c material is 
plotted in Fig. 4(a). As we can see, as Zn is introduced, 
T c reduces from 50K continuously down to 40K at 6% 
of Zn. The slow reduction in T c may suggest that the 
superconductivity saturates at large Zn doping. It will 
be interesting to confirm this by doping high Zn concen- 
tration under high pressure, which remains a challenge 
in material preparation. 

In addition to the critical temperature, another impor- 
tant feature during the transition from s± to s ++ is the 
change of low energy density of states (DOS). As shown 
in Fig. EJa), the gap amplitude reduces accompanying 
with the increase of low energy DOS when the system 



approaches the transition point tii 



0.04 from clean 



limit. And if one further increases the impurity concen- 
trations, the low energy DOS will be suppressed again 
due to the reopening of the gap. The non monotonic be- 
havior of DOS with impurities concentration should be 
able to be observed by integrated photoemission spec- 
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FIG. 4: (Color online) (a) T c versus Zn-doping concentration 
for SmFei-^Zn^AsOo.gFo.i- (b) Lattice constants a and c as 
functions of Zn-doping content in SmFei-^ZnzAsOo.gFo.i. 



(b) 400 




FIG. 5: (Color online) (a) DOS for various impurity concen- 
trations. n lmp = 0,0.015,0.025,0.035,0.045,0.065 and 0.15 
from bottom to top. (b) The low temperature specific heat 
at various impurity concentrations. 



troscopy. 

The change of DOS with impurity concentration may 
also be observed by the experiments which can measure 
the low energy DOS, for example the specific heat. In 
the superconducting state, the electron specific heat can 
be calculated with 



C(T) = ^ jT EN(E)f(E)dE, 



(11) 



where N(E) is the DOS, and f{E) is the Fermi distribu- 
tion function. In the low temperature regime, the tem- 
perature dependence of the superconducting order pa- 
rameter is very weak and can be neglected. So we use 
the zero-temperature DOS to calculate the specific heat 
with Eq. (fTTj) . In the calculation, we use A co h = 0.18ii = 
6meV as the energy scale, and the result is depicted in 
Fig. [5]Jb). We find that the electron specific heat below 
1QK is small in the clean limit, and shows a significant 



increase with approaching the transition point by increas- 
ing impurity concentrations. When the system is stabi- 
lized in the s ++ state, C(T) drops to a low value again. 
And the absolute value shown in Fig. [5jb) is in the same 
order or even lar ger than the experimental measurement 
of 1111 material [221. So it should be able to be observed 
in experiments. 

We note the recent work of Efremov et al. [23j, who 
applied T-matrix method to study the non-magnetic ef- 
fect on FeSC. Our microscopic theory shares some simi- 
larities with their. In their phenomenological theory the 
impurity-doping behavior is found to be associated with 
the averaged pairing coupling strength. In our theory, 
the decisive role of on-site pairing on the impurity effect 
is identified. 



V. SUMMARY 

In summary, we have studied the disorder-induced 
pair-breaking effect on the Fe-based superconductors us- 
ing a model that incorporates both the on-site and NNN 
pairings. We show that the Zn-impurity largely sup- 
presses NNN pairing. Its effect to the superconductivity 
depends strongly on the on-site pairing coupling strength 
go. The superconductivity can be robust, or evolves a 
transition from s± to s++, or is strongly suppressed in 
the presence of the disorder. Our theory qualitatively ex- 
plains different reductions of T c in various iron pnictide 
superconductors observed in the experiments on the Zn- 
impurity effect. We also predict the possible Zn-impurity 
doping induced transition from s± to s ++ pairing states 
in certain samples. Furthermore, we have systematically 
prepared Sm-1111 samples with T c = 50K under the Zn- 
doping and show that the reduction of T c could be con- 
sistent with the scenario of a moderate on-site s-wavc 
pairing. It will be highly interesting and important to 
prepare systematic controlled samples with higher Zn- 
doping to experimentally confirm or falsify the theory. 
The reduction of the gap in DOS during the transition 
from s± to s++ can be observed by integrated photoemis- 
sion spectroscopy or specific heat experiments. Finally 
we remark that the explicit paring forms are unlikely to 
be universal in Fe-based superconductors, in contrast to 
the universal d-wave pairing in cuprates. 
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